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Nuclear parton distribution functions (NPDFs) are determined by global analyses of experimen- 
tal data on structure- function ratios F^/F^ and Drell-Yan cross-section ratios ct^y/^dy- The 
analyses are done in the leading order (LO) and next-to-leading order (NLO) of running coupling 
constant a s . Uncertainties of the NPDFs are estimated in both LO and NLO for finding possible 
NLO improvement. Valence-quark distributions are well determined, and antiquark distributions 
are also determined at x < 0.1. However, the antiquark distributions have large uncertainties at 
x > 0.2. Gluon modifications cannot be fixed at this stage. Although the advantage of the NLO 
analysis, in comparison with the LO one, is generally the sensitivity to the gluon distributions, 
gluon uncertainties are almost the same in the LO and NLO. It is because current scaling-violation 
data are not accurate enough to determine precise nuclear gluon distributions. Modifications of the 
PDFs in the deuteron are also discussed by including data on the proton-deuteron ratio F® / F% in 
the analysis. A code is provided for calculating the NPDFs and their uncertainties at given x and 
Q 2 in the LO and NLO. 
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I. INTRODUCTION 

Cross sections of high-energy nuclear reactions are ex- 
pressed in terms of nuclear parton distribution functions 
(NPDFs), so that precise NPDFs are essential for find- 
ing any new phenomena in the high-energy reactions. Re- 
cently, this topic is becoming important in heavy-ion col- 
lisions for investigating properties of quark-hadron mat- 
ters [l| and also in neutrino reactions for investigating 
neutrino-oscillation physics 0] • Determination of precise 
NPDFs is valuable for studying various phenomena in 
heavy- ion reactions such as color glass condensate Q, 
J/ip suppression Q, and parton-energy loss The 
NPDF studies should be also important for heavy-ion 
collisions at LHC (Large Hadron Collider) @ . 

In neutrino oscillation experiments, most data are 
taken at small Q 2 (<1 GeV 2 ). We could approach such 
a kinematical region from the high-energy deep inelastic 
one by using quark-hadron duality 0- However, there 
are still unresolved issues in neutrino deep inelastic scat- 
tering. For example, an anomalous sin 2 9 W value was 
reported in the neutrino-iron scattering by the NuTeV 
collaboration |8|. It could be related to a nuclear mod- 
ification difference between the parton distributions u v 
and d v [1, [T(| [Tl| because the iron target is used in 
the NuTeV measurements. There is also an issue that 
nuclear corrections are different from the ones expected 
from electron and muon scattering experiments accord- 
ing to recent NuTeV data 

In these high-energy nuclear reactions, nucleonic PDFs 
rather than the nuclear ones are often used in cal- 



culating cross sections by neglecting nuclear modifica- 
tions although it is well known that nuclear corrections 
could be as large as 20% in medium-size nuclei 
These nuclear modifications have been experimentally 
investigated mainly by the measurements of structure- 
function ratios F^ / F^ and Drell-Yan cross-section ra- 
tios o~D Y 1 1(7 dy ■ Physical mechanisms of the nuclear cor- 
rections are, for example, summarized in Ref. [l3l |. In 
the small- a; region, the NPDFs become smaller than the 
corresponding nucleonic ones, which is called shadowing. 
There are depletions at medium x, which is related to the 
nuclear binding mechanism and possibly to a nucleonic 
modification inside a nuclear medium [141 ]. At large x, 
the nucleon's Fermi motion gives rise to positive correc- 
tions. 

Because the PDFs are related to the nonperturba- 
tive aspect of quantum chromodynamics (QCD), theo- 
retical calculations have been done by lattice QCD or 
phenomenological models. However, such calculations 
are not accurate enough at this stage. One would like 
to have accurate NPDFs, which are obtained in a model- 
independent way, for calculating precise nuclear cross sec- 
tions. We should inevitably rely on experimental data for 
determining them. 

Studies of nucleonic PDFs have a long history with 
abundant experimental data in a wide kinematical re- 
gion [IK]. However, determination of NPDFs is still at a 
premature stage with the following reasons. First, avail- 
able experimental data are limited. The experiments of 
the Hadron-Electron Ring Accelerator (HERA) provided 
data for structure functions at small a; in a wide range of 
Q 2 ; however, such data do not exist for nuclei. Because of 



final-state interactions, hadron-production data may not 
be suitable for the NPDF determination, whereas they 
are used in the nucleonic analysis. Second, the analy- 
sis technique is not established. Parametrization stud- 
ies for the NPDFs started only recently. The NPDFs 
are expressed in terms of a number of parameters which 
are then determined by a \ 2 analysis of the nuclear 
data. However, it is not straightforward to find func- 
tional forms of mass-number (A) and Bjorken-a; depen- 
dencies in the NPDFs. Furthermore, higher-twist effects 
could be important in the small-Q 2 region. 

A useful parametrization was investigated in Ref. (l6l | 
by analyzing F% structure functions and Drell-Yan data; 
however, the distributions were obtained by simply as- 
signing appropriate parameter values by hand in the ver- 
sions of 1998 and 1999. The first \ 2 analysis was re- 
ported in Ref. II 71 . and then uncertainties of the NPDFs 
were obtained [18[. All of these analyses are done in 
the leading order (LO) of the running coupling constant 
a s . A next-to-leading-order (NLO) analysis was recently 
reported [lj|. The LO x 2 analysis with the uncertain- 
ties was also investigated in the 2007 version of Ref. 
16 1 . There are related studies on the nuclear shadowing 
2C , Hi], [22j and a global analysis of structure functions 

M- 

In this way, the parametrization studies have been de- 
veloped recently for the NPDFs, and they are not suffi- 
cient. Here, we extend our studies in Refs. [l7], EH by 
focusing on the following points: 

• NLO analysis with NPDF uncertainties together 
with a LO one, 

• roles of NLO terms on the NPDF determination by 
comparing LO and NLO results, 

• better determination of x and A dependence, 

• nuclear modifications in the deuteron by including 
F2 (deuteron) /F2 (proton) data, 

• flavor asymmetric antiquark distributions. 

This article is organized as follows. In Sec. HH our 
analysis method is described for determining the NPDFs. 
Analysis results are explained in Sec. IIIII Nuclear mod- 
ifications in the deuteron are discussed in Sec. IIVI The 
results are summarized in Sec. W\ 



II. 



ANALYSIS METHOD 



The optimum NPDFs are determined by analyzing ex- 
perimental data of the F2 structure functions and Drell- 
Yan cross sections for nuclear targets. Details of our 
analysis method are described in Refs. [I3,[3> so that 
only the outline is explained in the following. 



A. Parametrization 

The parton distribution functions are expressed by 
two variables x and Q 2 . The variable Q 2 is defined 
Q 2 = —q 2 with the virtual photon momentum q in the 
lepton scattering, and the scaling variable x is given by 
x = Q 2 /(2Mv) with the nucleon mass M and the en- 
ergy transfer v. The variables for the Drell-Yan process 
are momentum fractions, x\ and Xi for partons in the 
projectile and target, respectively, and Q 2 defined by the 



dimuon mass as Q 



m. 



In our analysis, the NPDFs 



are expressed in terms of corresponding nucleonic PDFs 
multiplied by weight functions: 



f A (x,Q 2 )=w i (x,A,Z)f i (x,Q 2 ). 



(1) 



The functions f A and /, indicate type-* NPDF and nu- 
cleonic PDF, respectively, and W4 is a weight function 
which indicates a nuclear modification for the type-i par- 
ton distribution. The function Wi generally depends on 
not only x and A but also the atomic number Z. It should 
be noted that this expression sacrifices the large- x (> 1) 
nuclear distributions. Finite distributions exist even at 
x > 1 in nuclei, whereas the distributions should vanish 
in the nucleon. 

Flavor symmetric antiquark distributions are assumed 
for u, d, and s in the previous analyses [13, EB|- From 
the violation of the Gottfried sum rule and Drell-Yan 
measurements, it is now well known that these antiquark 
distributions are different [24| . It is more natural to in- 
vestigate modifications from the flavor asymmetric anti- 
quark distributions in the nucleon. In this work, flavor 
asymmetric antiquark distributions are used in nuclei, 
and the distribution type i indicates u v , d v , U, d, s, and 
9- 

u A (x, Ql) = w Uv (x, A, Z) ^Ei^) + Nd *( x > 



d A {x,Q 2 ) = w dv (x,A,Z) 
u A (x, Ql) = Wq(x, A, Z) 



d A (x,Q£) = Wq(x,A,Z) 

s A (x, Ql) = Wq(x, A, Z) s(x, Ql), 
g A (x, Ql) = w g (x, A, Z) g(x, Ql), 



A 

Zd v {x,Ql) + Nu v {x,Q 2 ) 
A _ 

Zu(x,Q 2 ) + Nd(x,Q 2 ) 
A 

Zd(x,Q 2 ) + Nu(x,Q 2 ) 



A 



(2) 



whereas the relation 



= d — s is assumed in Refs. 



< mi 



.17, 18]. The number of flavor is three (four) at Q 
(Q 2 > ml) [lq |. The bottom and top quark distri- 
butions are neglected. The strange-quark distributions 
are assumed to be symmetric (s — s) although there 
are recent studies on possible asymmetry s ^ s [25[. 
The charm-quark distributions are created by Q 2 evo- 



lution effects [18|. As for the nucleonic PDFs in the 



LO and NLO, the MRST (Martin, Roberts, Stirling, and 
Thorne) parametrization of 1998 is used [26[ in this anal- 
ysis, where the charm-quark mass is m c =1.35 GeV and 



2 



scale parameters are A4=0.174 and 0.300 GeV for the LO 
and NLO, respectively. In our previous analysis [Hj], the 
MRST-2001 version was employed. Since the NLO gluon 
distribution is negative at x < 0.005 and Q 2 —1 GeV 2 
in the 2001 version, we use the 1998 parametrization in 
this work. Negative gluon distributions in nuclei could af- 
fect our analysis inappropriately in the shadowing region. 
Furthermore, other researchers may use our NPDFs at 
small x (< 0.005) with Q 2 ~1 GeV 2 for calculating cross 
sections, for example, in LHC experiments. We tested 
various PDFs of the nucleon, but overall results are not 
significantly changed. Since we are interested in obtain- 
ing the distributions at Q 2 >1 GeV 2 in comparison with 
other distributions and also our previous distributions, 
we decided to use the MRST distributions. In the analy- 
ses of the CTEQ (Coordinated Theoretical/Experimental 
Project on QCD Phenomenology and Tests of the Stan- 
dard Model) collaboration, the initial scale Q 2 =(1.3) 2 
GeV 2 is used. 

The nuclear modification is assumed to have the fol- 
lowing functional form: 



Wi(x,A,Z) = 1+1 



1 

A a 



cii + biX + C{X 2 + d, L x 3 



(1 - x)fii 

where a, a%, bi, Ci, di, and /3, are parameters. The pa- 
rameter <2j controls the shadowing part, bi, Ci, and di de- 
termine a minute functional form, and Pi is related to the 
Fermi-motion part at large x. The parameter /3j is fixed 
at 0i=O.l because it cannot be determined from a small 
number of data in the Fermi- motion part. As it will be 
shown in the result section, the antiquark and gluon mod- 
ifications cannot be determined from the present data at 
x > 0.3. If a large pi, for example pi=l, is taken, the 
antiquark and gluon distributions could become unreal- 
istically large at large x. In order to avoid such an issue, 
Pi =0.1 is used. The overall A dependence in Eq. ((3]) is 
taken a = 1/3 simply by considering nuclear volume and 



surface contributions 27J ■ There are three constraints for 



the parameters by the nuclear ch arg e Z , baryon number 
A, and momentum conservation [1?], [H, [28[ : 

Z = Jdxj[2 u A (x, Ql) - d A (x, Ql)] , 

A = Jdxj [u A {x, Q 2 ) + d A (x, Ql)] , (4) 

A = J dxAx[u A (x,Q 2 )+d A {x,Q 2 ) 

+ 2 {u A (x, Ql) + d A (x, Ql) + s A (x, Ql)} + g A (x, Q 2 )] . 

We selected three parameters, a Uv , a^ v , and a g , which 
are determined by these conditions. 

Following improvements are made from the previous 
analysis [181 ]. First, the parametrization of £ dependence 
is modified. The meaning of the parameters bi, Cj, and di 
is not obvious, so that it is difficult to limit the ranges of 
these parameters in the analysis. Here, we take x points 
(x^, Xq { ) of extreme values for the function a, + biX + 



CiX 2 + diX 3 as the parameters instead of bi and Cj. They 
are related with each other by 



3di XQ { x Qi , 



3d; 



( X 0i + X 0i ) 



(5) 



Then, the values of x^ and x^ become transparent at 
least for the valence-quark distributions. From the mea- 
surements of the ratios F A jF® ', where D indicates the 
deuteron, we have a rough idea that the extreme values 



should be 



' Or 



0.15 and 



'Of 



0.6. The Droll- Yan 



measurements indicate x ^ ~ 0.1 for the antiquark dis- 
tributions. The value x^g and the extreme values for the 
gluon distribution are not obvious. 

Second, the A dependence 1/A 1 ^ 3 is too simple. From 
gross nuclear properties, the leading A dependence could 
be described by 1 — 1/A 1 ^ 3 . In order to describe more 
details, the parameters d Vl a,, dq, and d g are taken to be 
A dependent: 



= n« 1- 



V = V 



1 



A'r 



Tj d v , Qqi dg. 



(6) 



Because the extreme values a;i and Xq~ are almost in- 
dependent of A according to the F2 data and also from 
the previous analyses [13, HH, they are assumed to be 
independent of A. In addition, the parameters d^ and 
Xq 9 are fixed in the antiquark and gluon distributions as 
follows. Because the gluon distributions cannot be well 
determined from the present data, the parameter c g is 
taken c g = as assumed in the previous analysis [l8j . It 
means that Xq 9 and x$ g are related by Xq 9 = —Xq s from 
Eq. ([5]). There are still six parameters for the antiquark 
distributions and they should be too many in comparison 
with four parameters for the valence-quark distributions. 
We decided to fix the parameter , which is sensitive to 
the gluon shadowing ratio to the antiquark one because 
of the momentum conservation. We found that the gluon 
shadowing cannot be well determined even in the NLO 
analysis; therefore, the value of d^ is taken so that the 
gluon shadowing is similar to the antiquark shadowing. 

After all, the following twelve parameters are used for 
expressing the nuclear modifications: 

valence quark: Xq v , x^ v ,d^\ d^ , 

antiquark: x£ g , Xq 8 , af\ af\ df\ (7) 
gluon: z+.dW.dW 

These parameters are determined by the following global 
analysis. 

B. Experimental data 

Most of used experimental data are explained in Ref. 
[TH - First, the data for F A /Ff are from European 
Muon Collaboration (EMC) ^HaHH, the SLAC (Stan- 
ford Linear Accelerator Center)-E49, E87, E139, and 
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E140 collaborations [3ji [3j£ |3J, [35j] , the Bologna- CERN- 
Dubna-Munich-Saclay (BCDMS) collaboration [3(1 \3T\. 
the New Muon Collaboration (NMC) [11], the Fermilab 
(Fermi National Accelerator Laboratory)-E665 collabo- 
ration [H, S3, and the HERMES Second, the ra- 
tios F Al /F A2 (A 2 ^ D) are from the NMC [H,SiS3- 
Third, the Drell-Yan data are from the Fermilab-E772 
and E866/NuSea [IH] collaborations. 

Additional data to the HKN04 (Hirai, Kumano, Nagai 
in 2004) analysis are the ones for the deuteron-proton 
ratio F 2 /F^. These deuteron data are added because 
precise nuclear modifications are needed for the deuteron 
which is, for example, used in heavy-ion experiments at 
the Relativistic Heavy Ion Collider (RHIC) [J. The 
F% /Ft data are taken from the measurements by the 
EMC tH, the BCDMS [47[, the Fermilab-E665 [H, and 
the NMC [49[ . These data are used for extracting infor- 
mation on the flavor asymmetric antiquark distributions 
(u 7^ d) in the nucleon [24|. Therefore, the flavor asym- 
metric antiquark distributions in Eq. ^ are essential for 
a successful fit and for extracting information on modifi- 
cations in the deuteron. 



of the NPDFs and coefficient functions: 



+ C g {x,a s ) (£>g A (x,Q 2 ) 



(9) 



where Cn (x, a s ) and C g {x,a s ) are the coefficient func- 
tions [5lf . and is a quark charge. The symbol ® de- 
notes the convolution integral: 



(10) 



The proton-nucleus Drell-Yan cross section is given by 
the summation of qq annihilation and Compton processes 



da pA 



pA 



d °qq 



da pA 



dQ 2 dxF dQ 2 dxF dQ 2 dxF 



(11) 



d<Tgt 



Aira 2 



Because the DGLAP (Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi) evolution can be applied only in the per- 
turbative QCD region, the data with small Q 2 values 
cannot be used in the analysis. However, the data in dQ 2 dxF 9Q 2 s 
a relatively small-Q 2 region (Q 2 =l— 3 GeV 2 ) are valu- 
able for determining antiquark distributions at small x 
(x = 0.001 — 0.01). As a compromise of these conflicting 
conditions, only the data with Q 2 > 1 GeV 2 are used in 
the analysis. However, one should note that the data in 
the range, Q 2 =l— 3 GeV 2 , may contain significant con- 
tributions of higher- twist effects which are not considered 
in our leading- twist analysis. 



They are expressed in terms of the PDFs and subprocess 
cross sections: 



dyi / dy 2 



dQ 2 dxF dQ 2 dxF 



X [qi(yi,Q 2 W(y 2 ,Q 2 ) + qi(Vi,Q 2 )q?(V2,Q 2 )], (12) 



da pA _ Ana 2 
dQ 2 dx F ~ 9Q 2 s 



E 



( : / dyi dy 2 



da, 



dQ 2 dxp 
dQ 2 dxF 



(yi,Q 2 M(y2,Q 2 ) + qf(y 2 ,Q 2 )} 



h(yi,Q 2 ) + q l (yi,Q 2 )]g A (y2,Q 2 ) 



(13) 



C. x analysis 



The parameters are determined by fitting experimental 
data for the ratios of the structure functions F A and 
Drell-Yan cross sections. The total \ 2 



( odata 

E ( ' 



R 



thco\2 



^data^2 



(8) 



is minimized to obtain the optimum parameters. The 
ratio Rj ata indicates experimental data for F A /F A and 

CT w/ CT £>Y' an d -R} hco is a theoretical ratio calculated by 
the parametrized NPDFs. The initial scale Qo i s taken 
Qo=l GeV 2 , and the distributions in Eq. are evolved 
to experimental Q 2 points to calculate the x 2 by the 
DGLAP evolution equations All the calculations 

are done in the LO or NLO, and the modified minimal 
subtraction (MS) scheme is used in the NLO analysis. 
The structure function F A (x, Q 2 ) is expressed in terms 



The cross sections da qq -(n} and da qq -m indicate subprocess 
cross sections for qq annihilation processes in the LO and 
NLO, respectively. The da qg indicates the cross section 
for qg and qg processes. The NLO expressions of these 
cross sections are, for example, found in Ref. [52]. Effects 
of possible parton-energy loss in the Drell-Yan process 
[53| are neglected in this analysis. 

Using these expressions for the structure functions F A 
and Drell-Yan cross sections, we calculate the theoretical 
ratios Rf c ° in Eq. ©. The total X 2 is minimized by 
the CERN program library MINUIT. From this analysis, 
an error matrix which is the inverse of a Hessian matrix, 
is obtained. NPDF uncertainties are estimated by using 
the Hessian matrix as 

(14) 

where is the Hessian matrix, is a parameter, and 
£ indicates the optimum parameter set. The A% 2 value 
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TABLE I: Parameters obtained by the LO and NLO analyses. The parameters a Uv , a,d v , and a g are fixed by the three 
conservation conditions in Eq. @. Since they depend on nuclear species, they are explained separately in Appendix IA1 



Distribution 


4 


x 




d< 2 > 


a* 1 ' 


a< 2 > 


(LO) 














u A d A 


0.675 ± 0.005 


0.208 ± 0.004 


5.089 ± 0.315 


0.876 ± 0.143 


(Appendix) 


(Appendix) 




0.120 ± 0.007 


0.387 ± 0.027 


60 (fixed) 


0.260 ± 0.043 


-0.394 ± 0.027 


0.318 ± 0.038 


9 A 


8.35 ± 9.09 


—Xq (fixed) 


-799 ± 1561 


(0.194 ± 0.341) xlO" 5 


(Appendix) 


(Appendix) 


(NLO) 














A ,A 


0.629 ± 0.004 


0.178 ± 0.004 


6.157 ± 0.254 


1.063 ± 0.143 


(Appendix) 


(Appendix) 




0.239 ± 0.009 


0.119 ± 0.004 


150 (fixed) 


0.212 ± 0.019 


-0.540 ± 0.024 


0.295 ± 0.024 


9 A 


7.24 ± 3.97 


— Xq (fixed) 


-26.5 ± 27.8 


(0.109 ± 0.114) xlO -3 


(Appendix) 


(Appendix) 



determines the confidence region, and it is calculated so 
that the confidence level P becomes the one-a-error range 
(P = 0.6826) for a given number of parameters (N) by 
assuming the normal distribution in the multiparameter 
space. In the analysis with the twelve parameters, it is 
Ax 2 = 13.7. This Hessian method has been used for 
estimating polarized PDFs and fragmentation functions 
HH as well as nuclear PDFs in our previous version 
18fl • The details of the uncertainty analysis are discussed 
in Refs. [l8j, [54j as well as in the nucleonic PDF articles 
M. 



III. 



RESULTS 



A. Comparison with data 

Determined parameters are listed in Table Q] for both 
LO and NLO. Three parameters are fixed by the con- 
straints from baryon-number, charge, and momentum 
conservations in Eq. ([5]), and we chose a Uv , dd v , and a g 
for these parameters. The values of the obtained param- 
eters and their errors are similar in the LO and NLO. 
However, the errors indicate that there are slight NLO 
improvements in comparison with the LO results. 

Each x 2 contribution is listed in Table [TTJ The values 
suggest that medium and large nuclei should be well ex- 
plained by the current LO and NLO parametrizations. 
However, small nuclei are not so well reproduced. The 
LO fit (x 2 /d.o.f.=1.35) is better than the previous anal- 
ysis with x 2 /d.o.f.=1.58 [HI], which is partly due to the 
introduction of new parameters for the A dependence in 
Eq. ([6]). If the x 2 values of the LO analysis are com- 
pared with the ones in Table HI of Ref. [TH, we find 
that much improvements are obtained for the nuclear ra- 
tios, F 2 Li /F 2 D , F 2 C /F 2 D , F 2 e /F 2 , if* /if, Fp/F 2 c , 



2 > ± 2 
i pW i pBe 

and (JdyIv 

pFe i pD 



a 



DY 



1° 



DY ' 

and a 



rpSn I pC i^Pb I pC pCa / -phi „vCa i pD pFe i pBe 
/ 2 ' 2 12 ' 2 / 2 ' ® DY I ® DY i ° DY I a DY > 

whereas the fit becomes worse for 

/M " I dy/ (T d I y- The X 2 /d.o.f. is further re- 

duced in the current NLO analysis. 

According to Table QTJ the NLO results with 
% 2 /d.o.f.=1.21 reproduce the data better than the LO 
ones with x 2 /d.o.f.=1.35, especially in the following data 



TABLE II: Each x contribution. 



Nucleus 


Reference 


# of data 


2 

X 


2 

X 








ex r\\ 
(LO) 


( TVTT C\\ 

\r\YAj ) 


JJ/P 


\Aft A.7 AS 4Ql 


9Q0 


^7s a 
O ( o.o 


ozz .o 


4tt„ /ta 


35, 38| 


OO 


fin o 


01. o 


Lil/ D 


I38| 


1 7 
1 / 


oD.y 


oD.4 


D n /TA 

JoG/ L) 


|35| 


1 7 


oy.o 


Oo.U 


Kj/ D 


Ton on oo ac\\ 
Izy, oU, oo, oo, 4UJ 


4o 


lUo.o 


70 1-7 
to. ( 


IN / U 


[Of: /111 

1-50. 41] 


10Z 


loD.o 


1917 
Izl. ( 


a i /n 


loo, oo] 


"-t A 

oo 


^o.O 


AA Q 


Ca/D 


1 OU , OO, OO, 


33 


43.3 


34.1 




|OZ, O^, Oo, O / ] 


A 7 
'J / 


ins n 


Q7 A 


Cu/D 


[311 


19 


12.1 


13.2 


Kr/D 


[41] 


144 


115.1 


115.9 


Ag/D 


[35] 


7 


12.5 


9.1 


Sn/D 


[29] 


8 


13.3 


14.1 


Xe/D 


[39] 


5 


2.2 


2.3 


Au/D 


[34. 35] 


19 


55.6 


32.3 


Pb/D 


[40] 


5 


5.7 


4.5 


F 2 A /F 2 " total 




606 


792.4 


709.3 


Be/C 


[42] 


15 


12.6 


11.9 


Al/C 


[42] 


15 


5.0 


5.1 


Ca/C 


[38, 42] 


39 


29.9 


29.1 


Fe/C 


[42] 


15 


8.0 


8.3 


Sn/C 


[43] 


146 


204.0 


172.1 


Pb/C 


[42] 


15 


15.7 


12.2 


C/Li 


[38] 


24 


67.4 


64.9 


Ca/Li 


[38] 


24 


69.0 


65.3 


F Al /F A2 total 




293 


411.6 


369.0 


C/D 


[44] 


9 


9.3 


8.1 


Ca/D 


[44] 


9 


5.8 


13.8 


Fe/D 


[44] 


9 


12.6 


17.9 


W/D 


[44] 


9 


27.8 


29.6 


Fe/Be 


[45] 
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in this analysis to the data set of the previous version [ll| , 
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FIG. 1: (Color online) Comparison with experimental ratios 
R = Ff/Ff and F 2 D /ff . The rational differences between 
experimental and theoretical values {{R c * p ~ R thco )/ R thco ] are 
shown. The NLO parametrization is used for the theoretical 
calculations at the Q 2 points of the experimental data. The- 
oretical uncertainties in the NLO are shown at Q 2 =10 GeV 2 
by the shaded areas. 
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FIG. 2: (Color online) Comparison with experimental data of 
R = F^/Ff' Li . The ratios (i? cxp - R thco )/R thco are shown. 
The theoretical ratios and their uncertainties are calculated 
in the NLO. The notations are the same as Fig. [T] 



and they should provide a valuable constraint on PDF 
modifications in the deuteron. Because the F 2 jF 2 data 
are sensitive to u/d asymmetry |24j |. flavor asymmetric 
antiquark distributions should be used in our analysis. 
If the flavor symmetric distributions are used as initial 
ones, the fit produces a significantly larger % 2 . 

The fit results of the NLO are compared with the 
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FIG. 3: (Color online) Comparison with Drell-Yan data of 
R = ff£y/<r£y. The ratios (i? cxp - R thco )/R thco are shown. 
The theoretical ratios and their uncertainties are calculated 
in the NLO. The theoretical ratios are calculated at the Q 2 
points of the experimental data. The uncertainties are esti- 
mated at Q 2 =20 and 50 GeV 2 for the the ^r>Y/ a i^ type and 
(y V jjY I (y V jjY one > respectively. 



used data in Figs, [fl [2 and [3] for the ratios F 2 A /F 2 D , 
F 2 x jF 2 2 , and cr^/o^?, respectively. The rational 
differences between experimental and theoretical values 
(R cxp ~R thco )/R thco , where RisR = F^/Ff, F^/F^ 2 



or a 



pAi I pA 



DY 



1 1 a_ 



DY ' 



are shown. For the theoretical values, the 
NLO results are used and they are calculated at the ex- 
perimental Q 2 points. The uncertainty bands are also 
shown in the NLO, and they are calculated at Q 2 =10 
GeV 2 for the structure function F2 and at Q 2 =20 or 
50 GeV 2 for the Drell-Yan processes. The scale Q 2 =10 
GeV 2 is taken because the average of all the F2 data is 
of the order of this value. The scale is Q 2 =50 GeV 2 for 
the Drell-Yan ratios of the o^y/^dy type, and the lower 
scale 20 GeV 2 is taken for the ratio of the I& V dy type 



because experimental Q 2 values are smaller. 

These figures indicate that the overall fit is successful 
in explaining the used data. We notice that the x 2 val- 
ues, 53.0, 64.9, and 29.6 in the NLO, are especially large 
for F 2 Be /F 2 D , F§ '/F 2 Ll , and a^/a^y in comparison with 
the numbers of their data, 17, 24, and 9, according to Ta- 
ble [TTJ These large y 2 values come from deviations from 
accurate E139, NMC, and E772 data; however, such devi- 
ations are not very significant in Figs. [TJ[21 and [3] There 
are general tendencies that medium- and large-size nuclei 
are well explained by our parametrization, whereas there 
are slight deviations for small nuclei. Because any sys- 
tematic deviations are not found from the experimental 
data, our analyses should be successful in determining 
the optimum nuclear PDFs. 

Next, actual data are compared with the LO and NLO 
theoretical ratios and their uncertainties for the calcium 
nucleus as an example in Fig. |4j In the upper figures, 
the theoretical curves and the uncertainties are calcu- 
lated at fixed Q 2 points, Q 2 =10 GeV 2 and 50 GeV 2 for 
the F2 and the Drell-Yan, respectively, whereas the ex- 
perimental data are taken at various Q 2 values. The 
rational differences (R cxp ~ R thco )/R thco are shown to- 
gether with the difference between the LO and NLO 
curves, i? thoo (LO)/i? thco (NLO) - 1, in the lower fig- 



6 



EMC - E136 
• NMC • E66S 






Q 2 = 10 GeV 2 





FIG. 4: (Color online) Theoretical results are compared with 
the data of the F 2 ratio F° a /F° and the Drell-Yan ra- 
tio o v £y 1 / aV DY- I n the upper figures, the theoretical curves 
and uncertainties are calculated at Q 2 =10 GeV 2 for the F2 
ratio and at Q 2 =50 GeV 2 for the Drell-Yan ratio. The 
dashed and solid curves indicate LO and NLO results, and 
the LO and NLO uncertainties are shown by the dark- 
and light-shaded bands, respectively. The lower figures in- 
dicate the ratios (J? cxp - R thco )/R thco where R indicates 
F§ a IF® or (?z)Y /^dy- Here, the theoretical ratios are cal- 
culated at the experimental Q 2 points. For comparison, 
the LO curves and their uncertainties are also shown by 
J R thoo (LO)/ J R thco (NLO) - f and <5i? thco (LO)/7? theo (NLO). 



ures. The comparison suggests that both LO and NLO 
parametrizations should be successful in explaining the x 
dependence of the calcium data. It is noteworthy that the 
NLO error band of the F2 ratio becomes slightly smaller 
in comparison with the LO one at small x; however, mag- 
nitudes of both uncertainties are similar in the region, 
x > 0.02. The NLO improvement is not clearly seen in 
the Drell-Yan ratio o V dy I° V dy m the ran g e of x > 0.04. 
There are discrepancies between the theoretical curves 
and the F^jF® data at x < 0.01; however, they are 
simply due to Q 2 differences. If the theoretical ratios 
are calculated at the same experimental Q 2 points, they 
agree as shown in the Ft ? a /F® part of Fig. [T] 

These LO and NLO results indicate that the available 
data are taken in the limited x range without small- x 
data, and they are not much sensitive to NLO correc- 
tions. This fact leads to a difficulty in determining nu- 
clear gluon distributions because the gluonic effects are 
typical NLO effects through the coefficient functions and 
in the Q 2 evolution equations. 



B. Q 2 dependence 

One of possible methods for determining the gluon dis- 
tribution in the nucleon is to investigate Q 2 dependence 
of the structure function F2 [57j ■ Because Q 2 dependent 
data exist in the jF^ ratios, it may be possible to 
find nuclear modifications of the gluon distribution. 

The LO and NLO parametrization results are com- 
pared with Q 2 dependent data for F 2 Kr / ' F 2 D and F^ n jF§ 
measured by the HERMES and NMC collaborations, re- 
spectively, in Figs. O and [51 The uncertainties are shown 
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FIG. 5: (Color online) Comparison with Q 2 -dependent data 
of Ff r /Ff by the HERMES collaboration. The dashed and 
solid curves indicate LO and NLO results, and the NLO un- 
certainties are shown by the shaded bands. 
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FIG. 6: (Color online) Comparison with Q 2 -dependent data 
of Fi n /Ff by the NMC. The notations are the same as Fig. 



by the shaded bands in the NLO. Because the LO and 
NLO uncertainties are similar except for the small- a; re- 
gion, the LO ones are not shown in these figures. They 
are compared later in this subsection. The results indi- 
cate that overall Q 2 dependencies are well explained by 
our parametrizations in both LO and NLO. The compar- 
ison suggests that the experimental data are not accurate 
enough to probe the details of the Q 2 dependence. Fur- 
thermore, Q 2 dependencies of the HERMES and NMC 
results are different. The HERMES ratio F^/F® tends 
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FIG. 7: (Color online) Q 2 dependence of the ratio F 2 Ca /F 2 D is 
compared in the LO and NLO at £=0.001, 0.01, 0.01, and 0.7. 
The dashed and solid curves indicate LO and NLO results, 
and LO and NLO uncertainties are shown by the dark- and 
light-shaded bands, respectively. 



to decrease with increasing Q 2 at x=0.035, 0.045, and 
0.055, whereas the NMC ratio F$ n /Ff increases with 
Q 2 at the same x points, although the nuclear species 
are different. This kind of difference together with in- 
accurate Q 2 -dependent measurements makes it difficult 
to extract precise nuclear gluon distributions within the 
leading-twist DGLAP approach. It is reflected in large 
uncertainties in the gluon distributions as it becomes ob- 
vious in Sec. IIII CI 

In our previous versions [It], EH, the experimental 
shadowing in F^/F^ is underestimated at small x 
(0.01 < x < 0.02) partly because of an assumption on a 
simple A dependence. As shown in Fig. [5J the shadowing 
is still slightly underestimated at x = 0.0125; however, 
the deviations are not as large as before. If the exper- 
imental errors and the NPDF uncertainties are consid- 
ered, our paramctrization is consistent with the data. 

The NLO uncertainties are compared with the LO ones 
in Fig. 0for the ratio F 2 Ca /F 2 D . The LO and NLO ratios 
and their uncertainties are shown at 2=0.001, 0.01, 0.1, 
and 0.7. The differences between both uncertainties are 
conspicuous at small x (=0.001 and 0.01); however, they 
are similar at larger x. The LO and NLO slopes are also 
different at small x. These results indicate that the NLO 
effects become important at small x (< 0.01), and the 
determination of the NPDFs is improved especially in 
this small- a; region. 

Because the NLO contributions are obvious only in 
the region, x < 0.01, it is very important to measure the 
Q 2 dependence to pin down the NLO effects such as the 
gluon distributions. The possibilities are measurements 
at future electron facilities such as eRHIC [58| and eLIC 

M. 



C. Parton distribution functions in nuclei 

Nuclear modifications of the PDFs are shown for all 
the analyzed nuclei and 16 at Q 2 =l GcV 2 in Fig. [5] 
It should be noted that the modifications of u v are the 
same as the ones of d v in isoscalar nuclei, but they are dif- 
ferent in other nuclei. The modifications increase as the 
nucleus becomes larger, and the dependence is controlled 
by the overall 1/A 1 / 3 factor and the A dependence in Eq. 
([6]). The extreme values (x^, x^) are assumed to be in- 
dependent of A in our current analysis as explained in 
Sec. Ill A[ so that they are the same in Fig. [5] Although 
the oxygen data are not used in our global analysis, its 
PDFs are shown in the figure because they are useful 
for an application to neutrino oscillation experiments 0] ■ 
Our code is supplied at the web site in Ref. [60( for cal- 
culating the NPDFs and their uncertainties at given x 
and Q 2 . 

As examples of medium and large nuclei, we take the 
calcium and lead and show their distributions and uncer- 
tainties at Q 2 =l GeV 2 in Fig. [5] Because the deuteron 
is a special nucleus and it needs detailed explanations, its 
results are separately discussed in Sec. IIVI The figure in- 
dicates that valence-quark distributions are determined 
well in the wide range, 0.001 < x < 1 because the uncer- 
tainties arc small. It is also interesting to find that the 
LO and NLO uncertainties are almost the same. There 
are following reasons for these results. The valence-quark 
modifications at x > 0.3 are determined by the accurate 
measurements of F2 modifications. The antishadowing 
part in the region, 0.1 < x < 0.2, is also determined 
by the F2 data because there is almost no nuclear mod- 
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FIG. 8: (Color online) Nuclear modifications Wi (i — u v , d v , 
q, and g) are shown in the NLO for all the analyzed nuclei 
and 16 O at Q 2 =l GeV 2 . As the mass number becomes larger 
in the order of D, 4 He, Li, and Pb, the curves deviate from 
the line of unity (wi = 1). 
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FIG. 9: (Color online) Nuclear modifications of the PDFs and 
their uncertainties are shown for the calcium and lead nuclei 
at Q 2 =l GeV 2 . The dashed and solid curves indicate LO and 
NLO results, and LO and NLO uncertainties are shown by 
the dark- and light-shaded bands, respectively. 
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FIG. 10: (Color online) Nuclear modifications are shown for 
the calcium at Q 2 =100 GeV 2 . 



ification in the antiquark distributions according to the 
Drell-Yan data. If the valence-quark distributions are ob- 
tained at x > 0.1, the small- a; behavior is automatically 
constrained by the baryon-number and charge conserva- 
tions in Eq. ((3]). Because of these strong constraints, 
the LO and NLO results are not much different. In the 
near future, the JLab (Thomas Jefferson National Ac- 
celerator Facility) measurements will provide data which 
could constraint the nuclear valence-quark distribution 
especially at medium and large x (6l| . Furthermore, fu- 
ture neutrino measurements such as the MINERVA ex- 
periment at Fermilab [13] and the one at a possible neu- 
trino factory [63| should provide important information 
at small x (x = 0.01 - 0.1). 

The antiquark distributions are also well determined 
except for the large-a; region, x > 0.2, because there 
is no accurate Drell-Yan data and the structure func- 
tions F^ are dominated by the valence-quark distribu- 
tions. The LO and NLO uncertainties are similar ex- 
cept for the large-a; region. Because of gluon contribu- 
tions in the NLO, the antiquark shadowing modifications 
are slightly different between the LO and NLO. Pos- 
sible large- a; Drell-Yan measurements such as J-PARC 
(Japan Proton Accelerator Research Complex) [64[ and 
Fermilab-E906 65] should be valuable for the nuclear an- 
tiquark distributions in the whole- a; range. In the similar 



energy region, there is also the GSI-FAIR (Gesellschaft 
fur Schwerionenforschung -Facility for Antiproton and 
Ion Research) project [661 ]. 

The gluon distributions contribute to the F2 and Drell- 
Yan ratios as higher-order effects. Therefore, they should 
be determined more accurately in the NLO analysis than 
the LO one. Such tendencies are found in Fig. [5] because 
the NLO uncertainties are smaller than the LO ones in 
both carbon and lead. However, these NLO improve- 
ments are not as clear as the cases of polarized PDF [54| 
and fragmentation functions [55| . It is because the Q 2 - 
dependent data are not accurate enough to probe such 
higher-order effects as discussed in Sec. IIIIBI In order 
to fix the gluon distributions, accurate measurements are 
needed for the scaling violation of F£ /F^' [H, [HJ . The 
gluon distributions should be also probed by production 
processes of such as charged hadrons [63], heavy flavor 
68 1, T (69|, low-mass dilepton [70| . and direct photon 
711 ]. There is a recent study on the gluon shadowing 
from the HERA diffraction data [72j . Nuclear gluon dis- 
tributions play an important role in discussing properties 
of quark-hadron matters in heavy-ion reactions, so that 
they need to be determined experimentally. 

We also show the nuclear modifications at Q 2 =100 
GeV 2 for the calcium in Fig. [Til] The nuclear modifi- 
cations are not very different from those at Q 2 —l GeV 2 
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FIG. 11: (Color online) Parton distribution functions in the 
calcium at Q 2 =l GeV 2 . "No modification" indicates, for ex- 
ample, the distribution (Zu v + Nd v )/A in the figure of u„. 
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FIG. 12: (Color online) Parton distribution functions in the 
lead at Q 2 =l GeV 2 . "No modification" indicates, for exam- 
ple, the distribution (Zu v + Nd v )/A in the figure of u A . 



for the valence-quark distributions. However, the shad- 
owing corrections become smaller in the antiquark and 
gluon distributions in comparison with the ones at Q 2 =l 
GeV 2 , and the modifications tend to increase at medium 
and large x. 

The distribution functions themselves and their uncer- 
tainties are shown in Figs. [TTIandfT^lfor the calcium and 
lead, respectively. Both LO and NLO distributions are 
shown. Here, the uncertainties from the nucleonic PDFs 
are not included in the uncertainty bands. The calcium 
is an isoscalar nucleus, so that u A and u A are equal to d A 
and d A . However, they are different in the lead nucleus 
because of the neutron excess. In order to see nuclear 
modification effects, we show the distributions without 
the nuclear modifications. For example, the distribution 
(Zu v + Nd v )/A is shown in the figure of u A by using the 
MRST distributions for u v and d v . Although the uncer- 
tainties are large in the antiquark and gluon distributions 
at medium and large x, they are not very conspicuous in 
Figs. fTTI and [T^l because the distributions themselves are 
small. 

The flavor asymmetric antiquark distributions are as- 
sumed in this analysis as they are defined in Eq. 
From this definition, it is obvious that u A and d A are 
equal in isoscalar nuclei such as carbon and calcium. In 




0.001 



FIG. 13: (Color online) The ratio of flavor asymmetric dis- 
tributions, (u A — d A )/(u A — d A ), is shown for the proton, 
lithium, aluminum, iron, and lead at Q 2 =l GeV 2 . In the 
isoscalar nuclei, the distributions vanish {u A — d A — 0). 



Fig. [T3l the ratio (u A — d A )/(u A — d A ) is shown for the 
proton (p), lithium (Li), aluminum (Al), iron (Fe), and 
lead (Pb) at Q 2 =l GeV 2 . Because the nuclear correc- 
tions are assumed to be equal for the antiquark distri- 
butions at Qq in Eq. (JSJ) , they are almost independent 
of nuclear species except for the isoscalar nuclei. It is 
interesting to investigate possible nuclear modifications 
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on the distribution u — d at the future facilities [H, [6|| . 

There are noticeable differences between our NLO 
analysis results and the ones in Ref. [HI especially in 
the strange-quark and gluon modifications. These dif- 
ferences come from various sources. First, the analyzed 
experimental data sets are slightly different. Second, the 
strange-quark distributions are created by the DGLAP 
evolution by assuming s(x) = at the initial Q 2 scale, 
and the charm distributions are neglected in Ref. 
These differences lead to the discrepancies of the gluon 
modifications. 

The determined NPDFs and their uncertainties can be 
calculated by using our code, which is supplied on our 
web site [60l ]. By providing a kinematical condition for 
x and Q 2 and also a nuclear species, one can calculate 
the NPDFs. It is explained in Appendix [Bl If one needs 
analytical expressions of the NPDFs at the initial scale 
Q 2 , one may read instructions in Appendix [A] 



IV. NUCLEAR MODIFICATIONS IN 
DEUTERON 

Nuclear densities are usually independent of the mass 
number, which indicates that the average nucleon sepa- 
ration is constant in nuclei. However, the deuteron is a 
special nucleus in the sense that its radius is about 4 fm, 
which is much larger than the average nucleon separation 
in ordinary nuclei (~ 2 fm). Because it is a dilute sys- 
tem, nuclear modifications are often neglected. In fact, 
corrections to nucleonic structure functions and PDFs 
are small, namely within a few percentages according to 
theoretical estimates [zl [zl [zl, [zl, [n| even if they are 
taken into account. 

Theoretical modifications in the deuteron depend 
much on models. The shadowing of Ref. [731 ] . which was 
used in the MRST analysis, is 1.1% at x=0.015 and Q 2 =4 
GeV 2 ; however, other model calculations are different, 
0.5-3.5% at x ~ 0.01 as shown in Tabic EH Further- 
more, the modifications at medium x could be as large 
as or larger than the small- a; shadowing (for example, 5% 
at x = 0.6 — 0.7 according to Ref. [zl|). 

Such deuteron corrections are becoming important 
recently although the magnitude itself may not be 
large. For example, precise nuclear modifications need 
to be taken into account for investigating quark-hadron 
matters in heavy-ion collisions such as deuteron-gold 



reactions in comparison with deuteron-deuteron ones 
at RHIC 1]. They are also valuable for discussing 
Gottfried-sum-rule violation and flavor-asymmetric an- 
tiquark distributions because deuteron targets are used 

In our previous versions on the NPDFs [13, [Hj]> the 
data of /F% are not included in the used data set. 
Obtained nuclear modifications tend to be large in com- 
parison with the theoretical model estimations. For ex- 
ample, the HKN04 analysis [HI, [6(| indicates about 6 
and 8% corrections in antiquark and gluon distributions, 
respectively, at small x (~ 0.01) with Q 2 =l GeV 2 and 
about 5% in valence-quark ones at medium x (~ 0.7). 
They are possibly overestimations in the sense that typ- 
ical theoretical models have corrections within the order 
of a few percentages. 

In order to obtain reasonable deuteron modifications 
from experimental data, we added F^ jF^ measurements 
into the data set in our \ 2 analysis. These deuteron data 
were already included in the analysis results in Sec. IIIII 
In addition to the analysis of Sec. IIIII two other analyses 
have been made by modifying Eq. (|3|): 

/ 1 \ a l + b i x + ax 2 + diX 3 
Wi {x, A, Z) = 1 + e d ^1 - — j (1 _ x)0z • 

An additional factor e<j is introduced. The analysis re- 
sults in Sec. IIIII correspond to the ed=l case. The other 
analyses have been made by taking ed=0 and by taking 
it as a free parameter. We call them analyses 1, 2, and 
3, respectively: 

• analysis 1: e<j=l, 

• analysis 2: e<2=0, 

• analysis 3: ed=free parameter. 

The x 2 values of these analyses are listed in Table IIVI 

The deuteron modifications are terminated in the anal- 
ysis 2 by taking ed — 0. Although such an assumption 
does not seem to make sense, we found a smaller % 2 value 
from the / F^ data than the one of the first analysis 
as shown in Table HVl There are three major reasons for 
this result. First, the deuteron modifications obtained 
by the overall A dependence in Eq. ([3]) are too large 
because the deuteron is a loosely bound system which is 
much different from other nuclei. Smaller modifications 



TABLE III: Theoretical model estimations for nuclear modi- 
fication of F2 in the deuteron at x ~ 0.01. 



Reference Q 2 (GeV 2 ) x Modification (%) 

[73] 4 0.015 1.1 

[74] 4 0.010 2.5 

[75J 4 0.010 2.0 

[76] 4 0.010 0.4-0.8 

[77] 4 0.010 1-2 (3.5) 



TABLE IV: Each X 2 contribution. 



Nucleus # of data \ 2 je d = 1) x (e d = 0) \ 2 (free) 



D/p 


290 


322.5 


282.6 


284.9 


FtlF? 


606 


709.3 


704.8 


703.8 




293 


369.0 


381.6 


375.3 


Drell-Yan 


52 


85.1 


84.5 


85.4 


Total 


1241 


1485.9 


1453.4 


1449.4 


(xVd-o.f.) 




(1.21) 


(1.18) 


(1.18) 
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FIG. 14: (Color online) Comparison with the NMC data of 
/F%. The solid, dotted, and dashed curves indicate the 
NLO results of the analyses 1 (ej = 1), 2 (e<j = 0), and 
3 (ed=free), respectively. The shaded bands indicate uncer- 
tainties of the analysis-3 curves. 



are expected from the large nucleon separation. Second, 
deuteron data are used for determining the PDFs in the 
"nucleon" [26| by considering nuclear shadowing modifi- 
cations of Ref . [73] . The modifications are calculated in 
a vector-meson-dominance mechanism and the shadow- 
ing in the deuteron is about 1% at x ~ 0.01 according to 
this model. If this shadowing is not a realistic correction, 
the nucleonic PDFs of the MRST should partially con- 
tain nuclear effects at small x. The medium- and large-x 
regions are not corrected, so that some deuteron effects 
could be also included in the nucleonic PDFs in these 
regions. However, the corrections are not experimentally 
obvious in such x regions as we find in the actual data of 
Fig. [T3J In this way, the nucleonic PDFs could contain 
some deuteron modification effects. Third, the nuclear ef- 
fects could be absorbed into the u/d asymmetry because 
it is determined partially by the / F^ data. These are 
the possible reasons why analysis 2 produces the smaller 
X 2 value. 

In the third analysis, the additional parameter is 
determined from the global analysis. As mentioned, 
the internucleon separation is exceptionally large in the 
deuteron. It leads to small nuclear corrections, which 
are much smaller than a smooth A-de pen dent functional 
form, as calculated in various models [73l l74l (7f| [76|, [z3] 
in Table IIII1 These models could be used for estimating 
an appropriate value for e<j. However, we try to deter- 
mine the nuclear PDFs without relying on specific theo- 
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FIG. 15: (Color online) Nuclear modifications of the PDFs 
and their uncertainties are shown for the deuteron at Q 2 =l 
GeV 2 . The solid and dashed curves are obtained by analyses 
1 and 3, respectively, and their uncertainties are shown by the 
shaded bands. 



rctical models. The modification parameter e<j is deter- 
mined from the experimental data, and our analysis 3 
indicates e<j = 0.304 ± 0.141. If the diffuse deuteron sys- 
tem is considered, the 70% reduction may make sense. 
Nonetheless, it should be noted that this factor may not 
reflect a realistic modification because it is likely that the 
deuteron effects are contained in the nucleonic PDFs. 

We show actual comparisons with the experimental 
data for Ff/F% by the NMC in Fig. QH Three global 
analysis results are shown, and the shaded areas indi- 
cate uncertainty bands in analysis 3. The figure suggests 
that all the analyses are successful in explaining the data. 
However, as indicated in the x 2 reductions in the analy- 
ses 2 and 3, it is clear that their curves are closer to the 
experimental data at small x such as cc=0.005 and 0.008. 
There is a tendency that deviations from the data be- 
come larger as the deuteron modifications are increased. 
All the analysis results are more or less similar in the 
medium- and large- x regions. All the curves of the anal- 
ysis 1 (ed — 1) and 2 (e^ = 0) are within the uncertainty 
bands of the analysis 3 (e<2=free), although the analysis- 
1 curves are at the edges of the error bands at small x. 
It means that all these analysis are consistent with each 
other and they explain the experimental data. 

Modifications of the PDFs are shown for the deuteron 
in Fig. [15] at Q 2 =l GeV 2 . The results of analyses 1 and 
3 are shown with their uncertainty estimation. There is 
no deuteron modification in analysis 2 as obvious from 
the definition in Eq. (|15[) . The uncertainty bands of 
analysis 3 shrink at the points, where the nuclear mod- 
ifications vanish (w D = 1), for example, in the figure of 
the valence-quark modification w Uv . This is caused by 
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the error from the parameter e^. Its contributions to the 
uncertainties are large and the derivative dw Uv /ded is 
proportional to , which vanishes 

at the same points as the function w® v (x) — 1 = in Eq. 
(fT5|) . It leads to the gourd-shaped uncertainty band in 
w Uv because other terms are small. Such an error shape 
does not appear in analysis 1 because the error term of 
ea does not exist. Here, the derivative dw Uv / 'dSp is also 
a cubic polynomial and vanishes at three x points. How- 
ever, dw Uv /dxy v and dw Uv / dx^ v are quadratic functions, 
which do not vanish at the same x points, and their con- 
tributions to the uncertainties are of the same order of 
the dfcf 1 and dffl terms. This is the reason why such a 
gourd-shaped function does not appear in the uncertain- 
ties of the = 1 analysis. 

The antiquark shadowing is about 2% at x = 0.01 and 
the valence-quark modification is about 1% in the anal- 
ysis 1 according to Fig. Q1)J We should note that the an- 
tiquark shadowing is reduced about 30% at Q 2 =4 GeV 2 
as shown in Fig. in comparing it with the theoretical 
values in Table Hill In analysis 3, the antiquark shadow- 
ing becomes 0.5%, which could be slightly smaller than 
the theoretical ones. Both corrections (2% and 0.5%) 
are slightly different from the assumed correction (1% 
at £=0.01 and Q 2 =4 GeV 2 ) in the MRST fit. The un- 
certainty bands depend on the initial functional form or 
assignment of the parameters in the \ 2 analysis. The 
uncertainties are generally larger in analysis 3, and the 
line of w D — 1 is within the bands. It suggests that pre- 
cise deuteron modifications cannot be determined at this 
stage. The modifications in Fig. 1 151 may not be seriously 
taken because the deuteron effects could be partially in- 
cluded in the nucleonic PDFs. It is difficult to judge what 
the realistic deuteron modifications are at this stage. Ac- 
tual modifications are possibly in between these analysis 
results. 

In order to obtain realistic modifications, the nucleonic 
PDFs should be determined by considering the deuteron 
modifications, for example, of our analysis results. Then, 
using new nucleonic PDFs, we redetermine our nuclear 
PDFs including the ones in the deuteron. Realistic 
deuteron modifications should be obtained by repeating 
this step. 



V. SUMMARY 

Nuclear PDFs have been determined by the global 
analyses of experimental data for the ratios of the F 2 
structure functions and Drell-Yan cross sections. The 
uncertainties of the determined NPDFs are estimated by 
the Hessian method. The first important point is that 
the uncertainties were obtained in both LO and NLO so 
that we can discuss the NLO improvement on the deter- 
mination. We found slight NLO improvements for the 
antiquark and gluon distributions at small x (~ 0.001); 
however, they are not significant at larger x. Accurate 
experimental measurements, especially on the Q 2 depen- 



dence at small x, should be useful for determining higher- 
order effects such as the nuclear gluon distributions. 

The valence-quark distributions are well determined. 
The antiquark distributions are also determined at x < 
0.1; however, they have large uncertainties at x > 0.2. 
The gluon modifications are not precisely determined in 
the whole x region. Future measurements are needed to 
determine accurate nuclear distributions. 

Nuclear modifications were discussed for the deuteron 
in comparison with the experimental data of /F^. 
However, it is difficult to find accurate modifications at 
this stage because deuteron effects could be partially con- 
tained in the nucleonic PDFs. An appropriate nucleonic 
PDF analysis is needed in addition to accurate measure- 
ments on the ratio F® /F%, 

Our NPDFs and their uncertainties can be calculated 
by using the codes in Ref. [60| . 
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APPENDIX A: PARAMETERS a Uv , a dv , AND a g 

The determined parameters arc listed in TableQ] There 
are other parameters, a Uv , a,d vl and a g , which are au- 
tomatically calculated from the tabulated values by the 
conservation conditions of nuclear charge, baryon num- 
ber, and momentum in Eq. Because they depend 
on nuclear species, namely on A and Z, their calculation 
method is explained in the following. The flavor asym- 
metric antiquark distributions in Eq. ([2]) are used in this 
analysis, whereas the flavor symmetric ones are used in 
the previous versions in Refs. 0, EH, so that relations 
are slightly different from the ones in Appendix of Ref. 

One can obtain values of these parameters for any 
nuclei by calculating the following integrals: 



h(A) 
h 
h 



dx 
dx 
dx 



H v (x,A) 



(1 


-Xfv 




1 


(1 


- xf* 




X 



u v {x), I 2 (A) 



H v (x, A) . 

dx — nr- d v (x), 



u v (x), 



HA) 



dxx 



(1-x) 
H v (x,A) 
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(1 - x y» 
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(1- 


-xfv 
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(1- 
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X 


(1- 


-xfv 




(x,A) 



(i-xY 



d v (x), 
d v (x), 



9(x) 



(Al) 
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where f) v = /3g = fl g = 0.1 and Hi(x, A) is 
Hi(x, A) — biX + CiX 2 + di{A)x 



(A2) 



Here, the parameters d v , a g , d g , and d g depend on A 
according to Eq. (J6]). It is noteworthy that the integrals 
ii, I2, and I7 depend on A. Using these integrals, one 
obtains parameter values by the conservations in Eq. ([4J: 

a u { A,Z) = - Zh ^ + \ A - Z \^ A K (A3) 
" v ' ; ZI 3 + (A- Z)h 



B. Numerical values of the NPDFs can be obtained for 
given x, Q 2 , and A. In particular, if one is interested in 
estimating the uncertainties of the NPDFs, this practical 
code needs to be used. 



APPENDIX B: CODE FOR CALCULATING 

NUCLEAR PARTON DISTRIBUTION 
FUNCTIONS AND THEIR UNCERTAINTIES 



a g {A,Z) = -\- 
+ a dii (A,Z) ■ 



ZI 2 (A) + (A-Z)h(A) 
Zh + (A- Z)h ' 



a Uv (A,Z)^jI 5 



(A4) 









H) 







MA) 



(A5) 



From these values together with the parameters in Ta- 
bleland the nucleonic PDFs at Q 2 =l GeV 2 of the MRST 
parametrization [2(|, one obtains the NPDFs at Q 2 =l 
GeV 2 for a given nucleus. If distributions are needed at 
different Q 2 , one needs to evolve the NPDFs by using 
own evolution code. If one does not have such an evo- 
lution code, one had better use the code in Appendix 



Codes for calculating nuclear PDFs and their uncer- 
tainties can be obtained from the web site [6(|. A gen- 
eral guid eline for usage is explained in Appendix B of 
Ref. [18[ , and the conditions are the same in the current 
version, HKN07 (Hirai, Kumano, Nagai in 2007). There- 
fore, the details should be found in Ref. [T||. The NPDFs 
can be calculated for nuclei at given x and Q 2 , for which 
the kinematical ranges should be in 10~ 9 < x < 1 and 1 
GeV 2 < Q 2 < 10 8 GeV 2 . 

Input parameters for running the code are explained 
in the beginning of the file npdf07.f. A sample program, 
sample. f, is supplied as an example for calculating the 
nuclear PDFs and uncertainties. The uncertainties of the 
NPDFs arc estimated by using the Hessian matrix and 
grid data for derivatives with respect to the parameters. 
The usage is explained in the sample program. 
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